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' Abstract. Extreme Scattering Events and pulsar secondary spectra have highlighted funda- 
, mental problems in our understanding of the dynamics of interstellar turbulence. We describe 
04 ■ some of these problems in detail and present the theory behind the technique of speckle imaging, 
' which offers a prospect of revealing fundamental properties of the turbulence. It also offers the 
prospect of resolving pulsar magnetospheres on ~ 10 nas scales. 
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1. Why bother about turbulence? 

Compilations of pulsar scattering observations reveal that, in some average sense, the 
power spectrum of electron density fluctuations in the ISM follows a power law on scales 
(~| . from ^ 10 m up to lO^-i^m, with a power-law index close to 11/3, the value expected 
04 from Kolmogorov turbulence (Armstrong, Rickett & Spangler 1995). In the Kolmogorov 
Q . view of turbulence, mechanical and magnetic energy from macroscopic processes, such 

• as stellar winds and explosions, cascades to smaller scales via a succession of self-similar 
C/3 , kinetic-energy conserving turbulent eddies. It is surprising that the predictions of Kol- 

mogorov theory, which applies to incompressible hydrodynamic turbulence, should re- 
semble in any way the structure of the interstellar medium. 
T— I ■ The ISM is not incompressible over a large range of scales of interest (Luo & Melrose 

^ ] 2006), and is intrinsically ma(7neio- hydrodynamic. Indeed, the magnetic field must play 
a strong role in mediating the interstellar turbulent cascade because the presence of 
diffractive scintillation in pulsars shows that turbulence persists on scales many orders 

• of magnitude below the coUisional mean free paths of electrons and protons (Lithwick 
[ & Goldreich 2001). Moreover, there is no consensus on the most basic physics of how 

I ■ energy is mediated between large and small scales in the highly magnetized ISM, nor is 
_ there agreement on the interrelation between the turbulent velocity, magnetic field and 
I ■ density fiuctuations (Elmegreen & Scalo 2004, §4.12-13). 
^ ] Several key observations over the last decade have added a number of complications to 

•j-^ • our view of interstellar turbulence. The most notable has been the realization that there 
[ exist local pockets of anomalously strong turbulence. This is evident in the intcrmittency 
^ ■ of scintillations of intra-day variable quasars (Lovell et al. 2008). It is also manifest in the 
" " ' anomalous turbulence associated with tiny (^10^^ m) ionized clouds that are prevalent 
throughout the ISM and that are responsible for Extreme Scattering Events (Fielder et 
al. 1987). The discovery of strong parabolic arcs in the secondary spectra of many pulsars 
has also revealed that in many instances the scattering is often highly localized along the 
line of sight, and that the turbulence itself appears to be highly anisotropic. Attempts 
to incorporate these anomalous properties into a physical framework of the ISM have 
proven problematic and controversial (Spangler & Vazquez-Semadeni 2007). 
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We elucidate the problems exposed by ESEs and pulsar secondary spectra in §2, while 
§3 describes the method of speckle imaging that is allowing us to directly "image" these 
scattering structures in the ISM. In §4 we present some preliminary limits on the nano- 
arcsecond structure of the pulsars whose radiation is being scattered by some of these 
anomalous scattering structures. The final section describes the prospects for solving 
some of the fundamental questions that relate to interstellar turbulence. 

2. ESEs, Parabolic Arcs and Anisotropic scattering 

Extreme Scattering Events (ESEs) are intensity excursions exhibited by some compact 
quasars and pulsars and lasting between 10 and 50 days (Fiedler et al. 1987; Romani et 
al. 1987). The symmetric nature of their lightcurves argues that they are due to the 
passage of 4-70xl0^°m sized cloud-like features across the lines of sight. The observed 
event rate of 0.013 source"^ year~^ means that the clouds are common, with an estimated 
volume density of one per ~ 10~^pc^ (Fiedler et al. 1994; Walker & Wardle 1998). 

The existence of ESEs is problematic because their observed optical properties im- 
ply, at face value, internal pressures that exceed the diffuse ISM's by three orders of 
magnitude, assuming temperatures comparable to the diffuse warm ISM (Spanglcr & 
Vazquez-Semadeni 2007). In a simple model in which a plasma overdensity refracts the 
radio waves sufficiently to reproduce the caustic peaks observed in ESE lightcurves, 
column densities of ^ 10^^ cm^^ arc required. The volume density depends on the elon- 
gation of the structure along the line of sight; for an elongation i] = IOO772 the density is 
~ 10^r]2^ cm~^ (Romani et al. 1987). However, other models may explain ESEs without 
recourse to such extreme properties; Pen & King (2012) have recently proposed that ESE 
may instead be interpreted in terms of undeniense sheets in the ISM. 

There is evidence from pulsar scattering that much of the turbulence in the ISM is 
highly localised and anisotropic. One of the principal means of gleaning this information 
from pulsar scattering measurements is via the secondary spectrum. The secondary spec- 
trum, A(t,uj), is the squared amplitude of the two dimensional Fourier transform of the 
dynamic spectriim of a pulsar's intcinsity scintillations, I{i>, t). In the secondary spectrum, 
the conjugate of observing frequency is the delay, while the Fourier conjugate of time is 
Doppler frequency. In the regime of strong scattering from a thin scattering screen, one 
can represent the received wavefield as the sum of wavefields from a set of stationary 
phase points (or speckles) on the surface of the scattering disk, u = ajC**^ ^ with each 
stationary phase point possessing an amplitude aj and phase = (j){xj) + {xj — /3r)^/2rp, 
where 0(xj) is the phase delay imposed by the scattering medium at the position, x^, of 
the stationary phase point, and r is the location of the telescope on the observer's plane. 
The distance to the scattering screen, Dg, and the distance to the pulsar, Dp, also effect 
the total phase delay via the Fresnel scale, rp = {(iDs/kf''^, where /3 = 1 - Dg/Dp. The 
speckles emanate from positions Gj = (x^ — (3y) / Dg on the scattering disk. 

The intensity scintillation pattern, I{v,t) = uu* , is the result of the interference of 
each speckle with every other speckle on the scattering disk, and the resulting secondary 
spectrum takes the form (e.g. Walker et al. 2004), 

A{t,u)) oc '^ajak[S{T-Tjk)S{u>-L0jk)+S{T + Tjk)6{uj + Wjk)], (2.1) 

^jk = \ {dj - Ok) ■ Veff. (2.2) 

The location of power in the secondary spectrum is predominately dictated by the 
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Figure 1. Left: the secondary spectrum corresponding to 300 speckle points distributed accord- 
ing to isotropic Kolmogorov turbulence. The inset panel shows the spatial speckle distribution. 
Right: the same as the left figure except that the distribution of points has been made anisotropic 
by reducing the 6y co-ordinates of the speckles by a factor of 5. 

positions of the speckletfj. For any given pair of speckles, j and fc, power appears at 
both the co-ordinates {Tjk,i-Ojk) and {—Tjk^—ujjk). This symmetry is a reflection of 
the fact that A{t, lo) is derived from the Fourier transform of a real quantity, namely 
I{v^t). The effective scintillation velocity, Vcfi also influences the Doppler frequency of 
the speckles; it is usually dominated by the pulsar velocity, Vp, but it may in prin- 
ciple also be affected by the peculiar velocity of the screen, Vg or of the Earth, v®: 
veff = /3- Viss = [(1 - /3)vp + /3ve - v,] . 

The prescription given by eqs. (|2.1l) - (|2.2p affords a geometric interpretation of the sec- 
ondary spectrum. Consider interference caused by a highly elongated speckle pattern, in 
which all the speckles make a constant angle, a, to the scintillation velocity. Interference 
between a speckle at location Oj with the bright "core" of the image at 0^ = traces 
out a parabola with the locus (tq, wq) = [DsO^ /2c(3, OjVes cos a/A). Allowing variation in 
9k 0, we obtain the locus of an inverted parabolae whose apex occurs at co-ordinates 
(ro,a;o). If, instead, the distribution of speckles is not highly elongated, variations in a 
causes power to lie interior to a bounding parabola, as shown in Figure [TJ 

Many pulsars exhibit strong, sharply defined parabolic arcs of the sort seen in the right 
of Figure [TJ This reveals two important qualities of the scattering. 1 . The scattering is 
highly anisotropic. If it were otherwise, power in the secondary spectrum would instead 
lie interior to the main parabola. 2. The scattering occurs in localized patches, on a 
scattering screen that is thin, a few percent of the total pulsar distance. Variation in 
and P would otherwise smear out the locations of points in the secondary spectrum. 

Unfortunately, eqs. (|2.ip - (l2.2p do not enable an unambiguous reconstruction of the 
speckle distribution because the angle a is not known for any given speckle. The delay 
only measures 9^ and the Doppler frequency only measures ^jfcff cos a/A. This degener- 
acy can be broken by instead measuring the scintillations in the interferometric visibility 
on intercontinental baselines, as we now describe. 



3. Speckle Imaging 

It is possible to form a complete image of the distribution of speckles on the scattering 
disk by forming the secondary spectrum of the interferometric visibility. As we show 
below, the distribution of power in the visibility secondary spectrum is very close to that 
observed in the intensity secondary spectrum, but the visibilities provide information on 
the astrometric phase shift associated with each pair of speckles. In essence, one is able 

f The terms in square brackets in eq. 1)2. 2[) make only a small contribution to the overall delay. 
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to image the disk by using the secondary spectrum to isolate the wavefield of each pair 
of interfering speckles, and then measure the phase of this isolated wavefield from the 
visibility secondary spcctnim to localise the speckle positions to extremely high precision. 

To see how this works, consider the visibility that a two-element interferometer would 
measure from a scattered pulsar. Quantities measured with the first and second receiving 
elements are labelled with subscripts 1 and 2 respectively. The elements arc placed at 
locations ri = R — Ar/2 and r2 = R + Ar/2 on the observer's plane, with R the mean 
position of the telescopes, and Ar their relative displacement. The measured visibility is, 

N 

V{i^,t) = uxul = ^aija2,fc[cos($ij - $2,fe) + «sin($ij - $2,fe)], (3.1) 

where <I>jfe = — $2.fe is the phase difference between the jth stationary phase point 
measured at station 1 and the fcth stationary phase point measured at station 2. In the 
limit in which the bandwidth and observing duration are large, the Fourier-transform of 
the visibility dynamic spectrum reduces to, 

1 ^ 

1^(t,w) = — ^aija2,fe|exp[i$Pj<5(T + T,fe)'5(w + Wjfe)}, (3.2) 

j.k 

1 a($ij - $2,fc) ^ ,a a \ ^ ^ A 

^^•'^ = 2^ at = x^^^'- - ^''^^ • + A^^' • 

^% = cf>j-cf>k + ^ {elj - ei^,) - ^{drj + e2,k) ■ Ar, (3.5) 

F F 

where we write 6ij = (xj — /3R)/Ds, Ok = (x^ — (311) /Ds- Unlike its single-dish 
counterpart, the visibility secondary spectrum is not symmetric under the operation 
{t,uj) — >■ (— r, —u). This is because wjk 7^ — Wfej and Tjk ^ —Tkj- However, if the visibil- 
ities are measured on a baseline small compared to the scale of the scintillation pattern 
(i.e. Ar ^ DgGj), the amplitude of the visibility secondary spectrum is effectively iden- 
tical to that observed in the intensity secondary spectrum. 

An important difference between the intensity and visibility secondary spectrum is 
that the latter contains an cxp[i$°j.] phase term that is no longer antisymmetric. This 
is because the phase of each speckle measured at two widely-separated telescopes differs 
slightly because of astrometric phase term. In terms of the formalism introduced here, 
we see that when Ar = 0, the term proportional to the {6j + 6k) • Ar destroys the odd 
symmetry = —^kj- If '^^^ adds the contribution from the j,k and k,j terms in the 
secondary spectrum, it is possible to isolate the astrometric phase of each pair of speckles 
and determine the projection of the position of each pair of speckles, 6j + Ok along the 
baseline Ar. Thus, measurements along two baselines are sufiicient to to determine the 
position of each pair of speckles. This is the basis of scintillation speckle imaging. 

This technique was first applied to PSR B0834-I-06 by Brisken et al. (2010). These 
327 MHz observations revealed a number of intriguing properties of the scattering medium: 

• The scattering disk was composed of two separate structures, separated by 9 AU. 

• Speckles along the primary scattering disk, which was 16 AU long, were distributed 
anisotropically, with the ratio of the major to minor axis of the disk being at least 27:1. 
The distribution of speckles along the long axis of the primary scattering disk does not 
resemble that expected of Kolmogorov turbulence. 
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• The secondary scattering disk contributed about 4% of the total power. It is tempt- 
ing to speculate on the origin of this feature. One possibility, given the strong scattering 
properties the feature must possess in order to scatter a substantial amount of off-axis 
power back into the line of sight, is that it may be a cloud of the sort that is impli- 
cated in Extreme Scattering Events. This hypothesis, however, is difficult to test because 
ESEs are characterised by their optical properties when viewed on-axis with respect to 
a background source, whereas this object was off-axis by > 20 mas. 

3.1. Magnetic field limits 

Given the highly anisotropic scattering inferred towards PSR B0834-I-06 and other pul- 
sars which display strong parabolic arcs, it seems clear that the magnetic field plays 
an important role in the turbulent dynamics. One means of probing the magnetic field 
is to search for small rotation measure (RM) fluctuations associated with the scatter- 
ing medium. If RM fluctuations are present, the left and right-hand circularly polarized 
components of the pulsar radiation will experience slightly different phase delays, and 
this will cause different scintillations in each sense of circular polarization (Macquart & 
Melrose 2000), which would be visible in the circular polarization secondary spectrum. 

Brisken et al. (2010) report that, for PSR B0834-I-06, no detectable scintillating circular 
polarization signal was detected at the 0.1% level, and this places a limit on the RM 
difference of less than 1.2 x lO^^radm^^ across AU scales on the scattering disk. 

4. ISS as probes of nano-arcsecond pulsar structure 

The fact that speckle images exhibit structure across baselines of > 10 AU implies 
stringent constraints on the size of the pulsar emission region. Interference between the 
primary and secondary scattering disks in the case of PSR B0834-I-06 means that the 
radiation from the pulsar is at least partially coherent on a baseline oi 6 = 1.3 x 10^^ m 
as viewed at the scattering screen. This translates to a physical scale at the pulsar of 
5{Ds — Dp) = 4700 km. If the pulsar radiation had a gaussian angular brightness profile, 
the HWHM of the brightness distribution would be ^ 850 km. 

Although the radiation must be at least partially coherent on 9AU baselines, it is 
difficult to determine the degree of coherence. This is because one does not know what 
fraction of the total power should be received from the secondary scattering disk if the 
pulsar radiation were 100% spatially coherent on this baseline. One can in principle 
determine this by comparing the power associated with the interference between pairs 
of speckles on the primary disk, Pi^i, and pairs of speckles on the secondary scattering 
disk, P2-2, with the power associated with interference between speckles on the primary 
disk with those on the secondary disk, Pi-2- The pulsar is resolved if the quantity, 
TZ = Pi^iP2^2 1 Pi-2 significantly exceeds one. Figure [2] shows that Pi-i, P2-2 and P1-2 
can all be measured from the secondary spectrum. However, a complication arises because 
these three power quantities can be difficult to measure in practice, and it is difficult to 
relate 7^ to a specific measurement of the pulsar angular size. 

5. The Future 

There are two obvious prospects for progress in this field at present. The first, involving 
interstellar holography, which is not discussed in this short paper, has been advanced 
by the efforts of Walker et al. (2005, 2008). Holography takes advantage of the highly 
redundant information provided in the dynamic spectrum about the speckle distribution: 
for N speckles there are N{N —l)/2 interfering pairs measured in the secondary spectrum. 
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Figure 2. A schematic of tiie secondary spectrum of PSR 0834+06 found by Brislcen et 
al. (2010). Tiie green points represent interference between speckles on the primary scatter- 
ing disk, cyan represents interference between speckle pairs on the primary disk with those on 
the secondary scattering disk, and the purple points represent interference between adjacent 
speckles on the secondary scattering disk. 

The concept of interstellar holography may be viewed as a deconvolution problem. The 
Fourier transform of the intensity I{y,t) = u{v,t)u* {y^t), is just the autoconvolution of 
the Fourier-transformed wavefield: J(r, w) = u(t,w) -ku*{T,ui). 

Recent holographic work performed by Ue-Li Pen and collaborators on PSR B0834-I-06 
data has claimed a measurement of the pulsar's reflex motion (Pen et al. in prep.). 
Holography was used to effectively descatter the pulsar radiation and permit extremely 
high S/N measurements of the pulsar's radiation. By performing this holography over a 
succession of bins in pulse phase, it has been possible to measure a small but significant 
phase shift associated with the pulsar reflex motion. 

A second prospect is related to measuring the motions of speckles in scattered pul- 
sar images using a succession of speckle images over a period of weeks to months. In 
most cases, the pulsar proper motion dominates the effective scintillation velocity, but 
ultra- high S/N astrometric imaging, of the sort performed on PSR B0834+06, offers 
the prospect of resolving motions of the individual speckle groups relative to the bulk 
motion. This therefore offers the prospect of relating the density fluctuations associated 
with the scattering disk (which are, to some limited extent, recoverable with holographic 
techniques) with the underlying turbulent velocity fluctuations. 
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